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VQ ■ Abstract. We propose a new correlator in one-dimensional quantum spin 

chains, the s— Emptiness Formation Probability (s— EFP). This is a natural 
generalization of the Emptiness Formation Probability (EFP), which is the 
probability that the first n spins of the chain are all aligned downwards. In 
the s— EFP we let the spins in question be separated by s sites. The usual EFP 
$_^ , corresponds to the special case when s = 1, and taking s > 1 allows us to quantify 

"^T ■ non-local correlations. We express the s— EFP for the anisotropic XY model in 

a transverse magnetic field, a system with both critical and non-critical regimes, 
in terms of a Toeplitz determinant. For the isotropic XY model we find that the 
magnetic field induces an interesting length scale. 
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1. Introduction 

The Emptiness Formation Probability (EFP) is defined for a one-dimensional spin-1/2 
chain to be the probability of formation of a ferromagnetic string, i.e. the probability 
^^ ' that the first n spins of the chain are all aligned downwards, 

' P{n) = Z-iTr{e-^/'=«'^i^}, (1) 

f^ ■ where H is the Hamiltonian of the system, Z = Trjc^^/'^'^-^} is the partition function 

^ ! and 
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Here a^ is the usual Pauli spin matrix and j labels the sites on the chain. At zero 
i-rt ' temperature the EFP is given by 

g: P(n) = (0|F|0), (3) 

O ' where |0) denotes the ground state of the system. This fundamental correlator has 

been the subject of a good deal of recent interest 012 El E El El- It usually cannot be 
. , ■ calculated exactly, but a number of results have been obtained regarding its asymptotic 

r> ' behaviour as n — > oo. 

jrt ' In this Letter we propose a new correlator in one-dimensional quantum 

spin chains: the s-Emptiness Formation Probability (s— EFP). This is a natural 
generalization of the usual EFP. What we consider is the probability of finding n 
equally spaced sites aligned downward, for a given spacing s. Specifically, the s— EFP 
is defined as 

^.(")=(oin^^io)- (4) 

J=0 
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Notice that we take the Pauh matrices at the first n sites whose labels are multiples 
of s. The EFP thus corresponds to the particular case when s = 1. Since Q is 
not restricted to adjacent spins the s— EFP contains non-local information that is not 
available through the EFP. For the anisotropic XY model we shall derive its asymptotic 
form as n — > oo using the theory of Toeplitz determinants. For simplicity we take the 
temperature to be zero, but an extension to finite temperature is straightforward. 

2. EFP and s-EFP for the XY model 

Abanov and Franchini have recently considered the EFP for the anisotropic XY model 
in a transverse magnetic field h 0. They used the theory of Toeplitz determinants 
and extensions of what is generally known as the Fisher-Hartwig conjecture [7| |S] to 
obtain the asymptotic (n — > cx)) behaviour of P{n). The XY Hamiltonian is 

N-l . s N~l 

and it is critical along the lines j — 0,\h\ < 1 and h — ±1, i.e. quantum correlations 
decay algebraically when the parameters are on these lines and exponentially when 
they lie away from it ^ ^j . The EFP can in this case be written as 

Pin) = \detSn\, (6) 

where Sn is an ri x n Toeplitz matrix whose elements are 

(Sn),,^^ r e'^^-''^'aie)de. (7) 



The function cr(0) is called the symbol of the matrix 5„, and is given by 

1 1 cos 9 — h + i'y sin 6 , . 

'^W = :^ + 7. , ■ (8) 

^ ^-y/(cos6l-/i)2+72sin2e' 

Abanov and Franchini obtained the asymptotic form of the determinant lO in all 
regions of the 7 — ft, phase diagram. Their results are as follows jB]: for 1 > /i 7^ — 1, 
Piri) ~ E'e^"'^; for h > \ the previous result is multiplied by [1 + (— 1)"A]; on the 
critical lines h = ±1, P{n) ^ Err^l^^\V + (-/i)"A/-yrJ]e~"'^; and on the critical line 
7 = 0, |/i| < 1, P{n) ^ En~'^^^e~"' ", i.e. the decay is Gaussian. Explicit formulae 
were obtained for E, A, a and f3 as functions of 7 and h. The latter, for example, is 
given by 

P = -^ J ^ d9\og\aie)\. (9) 

For this system the s— EFP corresponds simply to 

Psin) = |det5„(s)|, (10) 

where Snis) is obtained from S'„ by removing the rows and columns whose labels are 
not multiples of s. Its matrix elements are 

(Snis))^,^:^ re''^'-''^'cT{e)de. (11) 
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This is in Toeplitz form, but cr(0) is no longer the symbol. To determine the symbol 
we must find the function a^ that satisfies 

a(a)e-'"""da = / as{a)e~"'°'da. (12) 

Multiplying H12|) by e™* with < 9 < 2tt, summing over n and using the Poisson 
summation formula, we arrive at 

The value of CTs at the point 6 is therefore obtained as the average value of the original 
symbol a over the s vertices of a regular polygon. Hence the s— Emptiness Formation 
Probability of the XY model ^ can also be computed as a Toeplitz determinant. 

Given the parameters 7 and h, the numerical calculation of the average (|13|l is 
a trivial task. Exact analytical results, on the other hand, are only available for 
particular cases. One of these cases is the isotropic XY model, which we analyze in 
the next section. 

3. The isotropic XY model 

As a simple application, let us consider the isotropic XY model, for which 7 = 0. In 
the non-critical regime |/i| > 1 the ground state is always a ferromagnet, i.e. the spins 
are all aligned up if ft. > 1, and all aligned down if /i < —1. Therefore P{n) vanishes 
in the first case and is unity in the second. Since the chain is translation- invariant, 
Ps{n) is the same as P{n) for any choice of s. 

The critical regime |/i| < 1 is of course much more interesting. The symbol (|Hl 
reduces to 

^ ^ f 1 if - fc < < fc, 
<t{9) = 1 - ' (14) 

I otherwise, 

where k is related to the magnetic field by 

h — cosfc. (15) 

Applying Widom's theorem |S] to the resulting Toeplitz determinant, one concludes 
that as n ^ cx) the EFP behaves like P{n) ^ En^^'^e"' ", where E and a are explicitly 
determined constants. In order to obtain the s— EFP we must compute the average 
P3|l . which in this case can be done explicitly. The result is a piecewise constant and 
even function, with jumps at the critical points ±9* given by 

[0, t:) 3 9* = min{a, 27r — a}, a — sk mod 27r, (16) 



and values 



s 



sk 

2^ 



e/s otherwise. 



where the brackets [•] denote the integer part and 

e = sign{a — vr}. (18) 

If s < n/k we have CTs(O) = 1/s and CTs(7r) — 0. In this case we get simply 
Ps{n) = s^^P{n), so the s— EFP decays much faster than the EFP, combining the 
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original Gaussian decay with an additional exponential term that depends on the 
spacing. On the other hand, if s > ir/k then the function as{9) is never zero, and 
instead of applying Widom's theorem we should apply the Fisher-Hartwig conjecture 
\7\. In that case the decay is no longer Gaussian, but becomes an exponential with a 
power-law prefactor, 

P,{n) - £;n-«e-^". (19) 

The constant E can be obtained from the Fisher-Hartwig formula and the other 
quantities are 

f3{s) = --r ln(a,(0)) + {7T- e*) ln(a, W)] (20) 

IT 

and 

e=^ln^(r7,(^)M(0)). (21) 

Interestingly, the decay of Ps{n) therefore shows a clear transition at the 
magnetic-field-dependent value s = ir/k, changing from Gaussian to exponential. In 
particular, for large spacings we see that both iTs(O) and iTs(7r) tend to k/n. Hence in 
this limit there is no power-law prefactor and the exponent saturates: 

lim /3(s) ==ln(7r/fc). (22) 

s — *oo 

It is important to observe that expression H2U|) for the exponent (3{s) is continuous 
when we consider s as a real number, despite the discontinuities that appear in (|17f) . 
The function as{0) is discontinuous whenever sk = 2n7r, but at those points 9* vanishes 
and hence /3(s) remains unaffected. The same is true for the jumps in the function 
(Ts(7r), which occur at sk = (2n -I- l)7r (due to the variable e), because then we have 
0* — TT. The derivative P'{s), on the other hand, is discontinuous: at the special points 
sk — m: the function /3(s) has a local minimum with a cusp form. Remarkably, for 
a fixed k its value at the minimum is independent of n, and is actually equal to its 
s —* oo limit, 

Pimr/k) = ln(7r/fc). (23) 

In Fig.l we plot the function /3(s) for different values of the magnetic field 
h — cos k. Since, of course, only integer values of s may be realized in the actual 
chain, we pick values of the form k = ir/i with an integer £. We see that this leads to 
the appearance of the above mentioned new length scale: /3(s) attains its minimal value 
whenever s is a multiple of £. It is worth remarking that Toeplitz determinants have 
also been used to describe quantum entanglement ^2E1E|) and that a generalization 
of the usual entanglement geometry, similar to the one proposed here, involving spins 
that are s sites apart, exhibits behaviour like that of the s— EFP |14j . 

4. The line 7^ + /i^ = 1 

In the phase diagram of the XY model the line 7^-1-/1^ ~ 1, which Franchini and 
Abanov call F^, is special. Along it the ground state is completely disentangled [T3|. 
i.e. it is given by 

|0)r. = Y[[cos{^/2)\ t), + (-1)^" sin(^/2)| i),], (24) 
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Figure 1. The renorinalized exponent /3(s)/ log(7r/A:), defined in 1201 . for some 
values of the magnetic field h = cos k. We see that for k = n/t there exists a 
typical length scale for the spacing s: the exponent is minimal for s = nl, and 
its value is equal to its asymptotic value. Notice that the exponent li(p) is not 
defined for s < i because in this regime the decay is Gaussian. The point s = £ 
marks a crossover in the behaviour of Ps(n). 



where | t)j and | l)j denote spin- up and spin-down states respectively at the site j, 
and the parameter i? is such that 



cos^ = 



1-7 

1 + 7' 



(25) 



Let us consider, for simplicity of notation, only the case h > 0. Then the symbol (jHJ 
is given by 



1 / z — cos^ 

2 \ 1 — z cos u J 

where z = e^^ . The average (fT^ becomes 

1 — COS1? sin^ t9 zcos""^!? 



2 1 — z cos* •& 



(26) 



(27) 



1 — cos i9 



1 + z cos* d 



and its modulus is 

l(^rj,(e)| , ^ ,, 

1 — z cos* V 

As expected, the s— Emptiness Formation Probability in this case is independent 
of s, because the exponent /3 is given by 

1 — cos -d ^ 



(28) 



/3 = -log 
and we thus have 

Ps{n) = 



1 — cos iS 



sin^"(z9/2), 



(29) 



(30) 



which is in fact an exact result. 
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